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Abstract: We study the <i-wave holographic superconductors (the <i-wave model 



(N 

proposed in [arXiv:1003.2991[hep-th]]) immersed in constant external magnetic fields 
. by using the analytic matching method and numerical computation. In the probe 

limit, we calculate the spatially dependent condensate solution in the presence of 
the magnetism and find that the expression for the upper critical magnetic field 
satisfies the relation given in the Ginzburg-Landau theory. The result shows that 
the upper critical field gradually increases to its maximum value B c2 at absolute zero 
temperature T = 0, while vanishing at the critical temperature T = T c . Moving 
away from the probe limit, we investigate the effect of spacetime backreaction on the 
critical temperature and the upper critical magnetic field. The magnetic fields as well 
as the electric fields acting as gravitational sources reduce the critical temperature of 
the superconductor and actually result in a dyonic black hole solution to the leading 
order. We obtain the expression for the upper critical magnetic field up to 0(k 2 ) 
order. The analytic result is consistent with the numerical findings. 
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1. Introduction 

The gauge/gravity duality [[I], |2], |3j has been proved to be a powerful tool for studying 
the strongly coupled systems in field theory. This duality provides a well-established 
method for calculating correlation functions in a strongly interacting field theory 
using a dual classical gravity description. The mechanism of the high temperature 
superconductors has long been an unsolved mysteries in modern condensed matter 
physics and the difficulties lie in the strong coupling nature of the theory. Consider- 
ing the above facts, Gubser first suggested that by coupling the Abelian Higgs model 
to gravity with a negative cosmological constant, one can find solutions that sponta- 
neously break the Abelian gauge symmetry via a charged complex scalar condensate 
near the horizon of the black holeQ, [5]. Later, Hartnoll et al proposed a holographic 
model for s-wave superconductors by considering a neutral black hole with a charged 
scalar and the only Maxwell sector A = A t . They captured the essence in this limit 
and showed that the properties of a (2+l)-dimensional superconductor can indeed 
be reproduced . The gravitational model that dual to the c?-wave superconductors 
was proposed in [0, § where the complex scalar field for the s-wave model is replaced 
by a symmetric traceless tensor. 

The purpose of this paper is to explore the behavior of the upper critical magnetic 
field for (i-wave holographic superconductors. We will work in both the probe limit 
and away from the probe limit. The probe limit corresponds to the case the electric 



- 1 - 



charge q — > oo or the Newton constant approaches zero. Away from the probe 
limit at a lower temperature, backreaction on the spacetime is important because 
the black hole solution becomes hairy and Coulomb energy of the matter field near 
the black hole horizon becomes larger. The phase diagram thus might be modified. 
In a more recent paper [0], an analytical calculation on the critical temperature of 
the Gauss-Bonnet holographic superconductors with backreaction has been presented 
and confirmed the numerical results that backreaction makes condensation harder [ITU, 
II], [12|, U, H 0> EE HE HI The hairy black hole solution requires 



to go beyond the probe limit. In ||, it was suggested to take finite q by setting 
2k 2 = 1. Recently, the author in || proposed to keep 2k 2 finite with setting q = 1 
instead. We will follow the latter choice. In p4[ , it was found analytically that for 
s-wave holographic superconductors the presence of the magnetic field results in the 
depression in T c , while the upper value of the critical magnetic field performance 
is improved (see [^5|, |26|, ^7], 25 for related work of analytic study on holographic 
superconductors). There have been a lot of works focused on the magnetic field 
effects on the holographic model of condensate matter EO, 130, BTI, 152, RBI B3, 1551, BB 



37| , |38| , |39| , [40f . In this work, we will generalize the analytic discussion on the upper 
critical magnetic field for s-wave superconductors in the probe limit [41] to the <i-wave 
holographic superconductors. We will consider the spacetime backreaction which has 
not been discussed in the available studies of the critical magnetic field for ci-wave 
holographic superconductors. Furthermore we will carry out numerical computations 
to check the analytic results for the critical magnetic field in d-wave holographic 
superconductors. We will check whether the result found in [fM| is universal which 
can hold in the <i-wave holographic superconductors. 

The organization of the paper is as follows: we first study the condensation 
of the order parameter in the probe limit in section 2. In the weak field limit, 
the critical temperature and the order parameter operator will be calculated first. 
Then we continue the calculation to the strong field limit and obtain an analytical 
expression for the backreaction on the upper critical magnetic field. In section 3, we 
study the effect of spacetime backreaction on the critical temperature and the upper 
critical magnetic field. The presence of the magnetic field actually leads to a dyonic 
black hole solution and the critical temperature drops due to the backreaction of 
the magnetic field. In section 4, we show at qualitative level, the analytic study is 
comparable with the numerical computation. The conclusion will be presented in 
the last section. 



2. The probe limit 

In this section, we will study the condensate solution and the upper critical magnetic 
field for d— wave superconductors by using the matching method. Up to now, two 
types of d— wave superconductors have been proposed: One is constructed by using 
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a symmetric, traceless second-rank tensor field and a U(l) gauge field in the back- 
ground of the AdS black hole|Q. The other model for the d— wave order parameter 
is dual to a charged massive spin two field propagating in an asymptotically AdS 
geometry ||. We will take the first one as the example. 

2.1 The critical temperature 

We adopt the action for the (i-wave superconductor as follows 

S= ^ l rf V^{(#+^)+£m}, (2-1) 



c =-1 



(D^YD^B^ + m 2 B^B^ + \f„ v F» v 



(2.2) 



where B^ v is a symmetric traceless tensor, = + iA fl is the covariant derivative, 
q and m 2 are the charge and mass squared of B^ u , respectively. The <i-wave super- 
conductors can condensate on the x — y plane on the boundary with translational 
invariance, and the rotational symmetry is broken down to Z(2) with the condensate 
changing its sign under a tt/2 rotation on the x — y plane. Same as in 0, we use the 
spatial dependent ansatz for the (i-wave superconductors 

B^ = diagonal(0, 0, ip{z), -i/>(z)), A = </>(z)dt. (2.3) 

In the probe limit, the bulk gravitational theory is described by the AdS-Schwarzschild 
metric 

ds 2 = £ (^-f(r)dt 2 + ]T dxfj + ^dr 2 , (2.4) 
where the metric coefficient 

MI 2 r 3 

/(r) = l-q- = l-^±, (2.5) 



and I is the AdS radius and M is the mass of the black hole. The Hawking temper- 

3M 1 / 3 

4^F73 ■ 



ature of the black hole is T = Setting z = — , the metric can be rewritten in 



the form 

ds 2 = l ^ [~f(z)dt 2 + dx 2 + dy 2 } + -^dz 2 , (2.6) 
z 2 L J z z j(z) 

where 

f(z) = l-z\ a = r -± = -irT. (2.7) 
The equations of motion of the two field ip(z) and 4>(z) are given by 

w + if + -w + — * - %^ - %^ + 4W = °> ( 2 - 8 ) 

z z z z z z a z j 

<t>" - = °- ( 2 -9) 
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The equations of motion of the d— wave superconductors are very similar to the 
s— wave model and the matching method should be very efficient. Before solving the 
above equations, let us impose the boundary condition near the horizon and in the 
asymptotic AdS region, respectively: 

1) . On the horizon z = 1, the scalar potential must be vanishing = and ip should 
be regular. 

2) . In the asymptotic AdS region z — > 0, the solution of the scalar field behaves like 

*P~C A _z A - +C A+ z A + } (2.10) 

where A± = - 1± ^ l7 +^ vfli m The coefficients C*a_ represents as the source of the dual 
operator and C A+ correspond to the vacuum expectation values of the operator that 
couples to B^v at the boundary theory. We require m 2 l 2 > —4 (thus A + < 0) 
such that the C A+ term is a constant or vanishing on the boundary. Note that the 
Ca_-2 A ~ term does not impose a constraint on m 2 l 2 by requiring that the third term 
on the left hand side of Eq.( [2.8p to be smaller than the other two terms since we 
have imposed C A _ = 0. The condensate of the scalar operator O in the boundary 
field theory dual to the field B^ u is given by 

<o„>=hf _ C V |. (2.n; 



In what follows, we choose to set the mass of the B^ v to be m 2 l 2 = — ~, so that A 
and l 
form 



fiv ^ hi, u — 4 , du uiiau — 2 

and A_ = — |. The asymptotic value of the scalar potential at the boundary has the 



<j ) (z)=fi-qz J (2.12) 

where q = p/r+. Here \i is interpreted as the chemical potential and p as the charge 
density in the boundary theory. 

Now we are going to solve the equations of motion by using the analytic method 
developed in |^2]] . Expanding the two field if) and <fi near the horizon 

<j>{z) = 0(1) - 0'(1)(1 -z) + i^(l)(l - zf + ... (2.13) 

i>{z) = m - -z) + Inm - zf + ... (2.14) 

and noting that regularity at the horizon which gives 

23 

^'(1) = "12^(1), (2-15) 



the expression for ip"(l) and 0"(1) can be derived from fl2.8|) and ( p.9|) , respectively 

288 rv ' 18a 2 



^ , ( 1 ) = ^^(1)-t^^(1), (2.16) 



</'(!) = (2-17) 
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The approximate solutions for ip and (j) near the horizon can then be written as 

4>{z) = -0'(1)(1 -z)- ^V(l)(l - z)\ (2.18) 
= m + ^(l)(l - z) + U^- - ^1 - ,) 2 . (2.19) 



12 rv /v 7 2 V 288 18a 

Connecting the near horizon solutions ( 2.18 ) and ( |2.19| ) with the boundary solutions 
( |2.10D and ( |2.12| ) at the intermediate point z m = 1/2 smoothly, we find 

q b a 2 b 
2 ~ 2 + 6a 2 
2a 2 b 

3a 1 



^-^=0 + 7^' ( 2 - 20 ) 
= (2.21) 



+ v 2 ; V2304 144a 2 / ' v ; 

2 A+ V V 576 36a 2 J ' V ; 

where we have defined —<p'(l) = b and = a. From ( |2.2U| ) and ( |2.21| ), we obtain 

b 3 . . 

^ = 4 + 4?> ( 2 - 24 ) 



3 /?_.2 /i & 



2 V6 a V 1 -? <2 ' 25) 

In order to evaluate the expectation value of the operator < Os >= \/2CA + r+ 2 , we 
eliminate the a& 2 term from (2.21) and ( |2.23|) and obtain 

C A+ = -^LV2a. (2.26) 
+ a/42 



For non-vanishing a, we can eliminate Ca + to obtain 



By further using the relation q = a = f 71 "^, the expectation value of the operator 
< 03 > is given by 



, = ^™^l + yi-%, (2.28) 
where the critical temperature is defined as 

T c = 27r( ^ )1/4 * 0.058VP- (2.29) 

The expectation value give in Q2.28 ) shows us that the d— wave condensate is indeed 
the second order phase transition with the mean field critical exponent 1/2. 



- 5 - 



2.2 The upper critical magnetic field 



Now we consider the case when the d— wave superconductor is immersed in a strong 
external magnetic field. To the leading order, we consider the ansatz 



ip = ip(x,y,z), 4> = (f) (z) J A y = B c2 x. 
The equation of motion for the B^ field now becomes 



(2.30) 



2/ 



1 



2iA 



or 



a' 



z z z 



of 



« 2 / 



i/> = 0, 



(2.31) 



where the prime ' denotes the derivative with respect to z. By assuming ip 
e wy F{x, z;p), the above equation of motion can be changed into 



z z a z j z A z A 



F(x,z;p) 



a- 



-d 2 x + (p-B c2 x) 



F(x,z;p). 
(2.32) 

This equation can again be solved by separating F as F(x, z;p) = X n (x;p)R n (z). 
We then obtain the following eigen equations 



U 2 X 
-du + —)X n (x;p) = Y x n(x;p), 



R" + ( y- + - ) R' n 



m 2 2 | 4 2f | B c2 \ r , 
z 2 f a 2 f 2 z 2 zf a 2 f 



Rn(z), 



(2.33) 
(2.34) 



where U = y / 2B c2 (x — Eq. (|2.33| ) can be solved by the Hermite polynomials as 
follows 

X n (x;p) = e- u2 ^H n (x), (2.35) 

where A n = 2n + 1 is the corresponding eigenvalue and n = 0,1,2... denotes the 
Landau energy level. We will focus on n = case, which corresponds to the droplet 
solution. Now we solve the equation ( [2.34 ) in the strong field limit. Firstly, let us 
expand R(z) near the horizon 



Ro(z) = R (l) - R' (l)(l -z) + - z) 2 + ... 



Note that the regularity at the horizon gives 



f 2 MD - ||«o(i). 



On the other hand, near the AdS boundary z — > 0, it sets 

Ro(z) = C A+ z%. 



(2.36) 



(2.37) 



(2.38) 
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From ( j2.34|) , the second order coefficients of Ro(z) can be calculated as 

*KD = ^*(U + ^?«o(D + ^MD - ^ft(D (2.39) 
Finally, we find the approximate solution near the horizon 

r. , x ^ /,x /23 5 c2 \ „ /<w , 1 /1897 5 c2 2 47£ c2 

RJ Z ) = R (l) +(_ + _££ i? 1 1 - + - + —^7 + £ 

ov ; ov ; V 12 3 " J 2 V 288 18a 4 36a 2 

,/(1,J WA.iDil :) 2 (2.10) 



18a 2 

Now let us match the solutions ( 2.38|) and ( p. 40 ) at the intermediate point z m = 1/2. 
Requiring the solutions to be connected smoothly, we have 

1 3/2 _ /6409 E c2 2 6 2 955 c2 \ 
Ca+( 2 } " ^ + il4^~M4^ + 288^J i?o(1) ' (2 ' 41) 



2 +v 2 y V 576 36 « 4 36 « 2 72 « 2 , 
From the above equations, we find the solution for b 2 = |0'(1) 2 | 

, (// , 2| 569 „ B 2 31231 2 , A , 

|0'(1) Hn^ 2 + ^ + ^-« - ( 2 - 43 ) 

When the external magnetic field _B c2 is vanishing, the above result return to 02.271) , 
which is crucial for the expression of the critical temperature given in ( [2.29 ). By 
plugging |0'(i)l = -^f, a = ^j- and (|2.29| ) into Q2.43|) , we obtain the ansatz for B c2 




B c2 = ^T 2 ( x 1 105144^ + 218617 - 569^) (2.44) 

Figure [I] shows the temperature dependence of the upper critical magnetic field, 
which is consistent with the Ginzburg-Landau theory. We know that for type II su- 
perconductors, in the upper critical field the cores of the vortices are nearly touching 
and the flux contained in each core is only one quantum flux, the average magnetic 
field is then B c2 ~ n$°T) • ^ ms sense 5 we can determine the superconducting 
coherence length 

t(X) « (1 - h~ Xl \ ( 2 -45) 

■L C 

where the critical exponent —1/2 is what wanted for £(T) in the Ginzburg-Landau 
theory. 



3. Backreaction 

Away from the probe limit, the background spacetime becomes non- neutral and 
hairy. Both the electric field and the magnetic field could backreact on the back- 
ground geometry. In this section, we first consider the backreaction of the electric 
field and then we turn to discuss the backreaction of the external magnetic field. 
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Figure 1: (color online) The upper critical magnetic field changes with T/T c for d— wave 
superconductors. T c denotes the critical temperature without external magnetic field. 



3.1 The critical temperature with backreaction 

The hairy black hole solution is assumed to take the following metric ansatz 

Aw*^ 1 it* 

ds 2 = -f(r)e-^dt 2 + ^ + ^(dx 2 + dy 2 ), 

fir) 1 

The tt and rr components of the background Einstein equation yield 

,ip 2 (f) 2 e x 



(3.1) 



f, , / 3r , 2 

r l z 



e //2 i 2 / 2 i ft V 7 ° " . ;/2\ 

— + m ip + /(- 



./' 2 



X' + 2/? 2 r(^' 2 + 



"7^ - 



o. 



(3.2) 
(3.3) 



When the Hawking temperature is above a critical temperature T > T c , the solution 
is the well-known AdS-Reissner-Nordstrom black holes 



./ 



T 



1 f r 



K 2 P 2 



I 2 2r A 



K 2 p 2 
~2^' 



X = if) = 0, = p 



1 



(3.4) 



Near the critical temperature T ~ T c , the AdS-Reissner-Nordstrom solution becomes 
unstable against perturbation of the scalar field because the coupling of the scalar to 
gauge field induces an effective negative mass term for the scalar field. This negative 
mass term will drive the scalar field tachyonic as the temperature is lowered at fixed 
charge density. At the asymptotic AdS boundary (r — > oo), the scalar and the 
Maxwell fields behave as 



V> = — -a + — -x^—' 



P . 
p h ... 

r 



(3.5) 
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where \x and p are interpreted as the chemical potential and charge density of the 
dual field theory on the boundary. 

The equations of motion for field and corresponding Maxwell field are 

i> + (j - y - + (-jr " y " jy)*P = 0, (3.6) 

0" + ( 2 - + |)0' -^0 = 0- (3.7) 

Near the critical temperature, the order parameter is small-valued and one can con- 
sider it as an expansion parameter 

e =< O a+ > . (3.8) 

It is worth noting that given the structure of our equations of motion, only the even 
orders of e in the gauge field and gravitational field, and odd orders of e in the scalar 
field appear here. So, we expand the scalar field ip, the gauge field as series in e 

= 0o + e 2 02 + e 4 04 + -- (3.9) 
^ = e^ 1 + e 3 ^ 3 + e 5 ^5 + - (3.10) 

The background metric elements f(z) and x{ z ) can De expanded around the AdS- 
Reissner-Nordstrom solution 

f = fo + e 2 f 2 + e A f 4 + ... (3.11) 
X = e 2 X2 + e 4 x 4 + ... (3.12) 

The chemical potential \i should also be corrected order by order 

fj, = pio + e 2 ^2, (3.13) 

where 8/j,2 is positive. Therefore, near the phase transition, the order parameter as 
a function of the chemical potential, has the form 



V -fJ>o \ 

5fi2 J 



1/2 



(3.14) 



We can see that when \i approaches /io, the order parameter becomes zero and phase 
transition happens. So the critical value of /i is fx c — jiQ. The critical exponent 1/2 
is the universal result from the Ginzburg-Landau mean field theory. The equation 
for is solved at zeroth order by O — A*o(l — z ) an d this gives a relation p = Pqt + . 
Thus, to zeroth order the equation for / is solved as 

m = 4 ( 1 - z) (l + z + , 2 - . (3.15) 



zH 2 \ \ 2H 
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Now the horizon locates at z — 1. We will see that the critical temperature with 
spacetime backreaction can be determined by solving the equation of motion for ip 
to the first order. Using matching method at the horizon, we obtain 



2 2 



7o 2 (i) 



-2)Vi(l). 



At the asymptotic AdS boundary, we have 
We expand ipi m a Taylor series 



(3.16) 



(3.17) 



fa = Vi(l) - ^(1)(1 -z) + -Vi(l)(l - *) 5 



(3.18) 



The second order of ■?/>'/ can be obtained by the equation (|3.6| ) 



<(l) = -(5 + 



2 2 

rim 



2/^(1) 2/^(1)2 
Then the expression of t/>i can be rewritten as 



)V4(i)- 



V>i(z) = Vi(l)-(2 + 



2 2 

rim 



/o 2 (l) 



/6'(i) , ^' 2 



fn(l) rim 2 , . rim 2 
+ 3 + ( 5 + 777777T ~ TTTTTT^X 2 + 



2/6(1) 2/6(1)' 



/o 2 (l) 



/o(l) 2/6(l) 2 
^i(l)(l-^) 2 



(3.19) 



Since the field function and its first derivative should be connected smoothly, we can 
properly match the expressions of field function at a middle point 



7"^ 777 ^ I 

&C + = ^(1) - (2 + ^)^(l)(l - z m ) - - 



3 , fSO) , ^ 



/6(1) 2/^(1)' 



+ (5 + 



/.?(!) 



2 2 
rim 



)(2 + 



2 2 
rim , 



2/6(1) 2/6(i) 2 ^ / 2 (in 



^i(l)(l -z m f 



(3.20) 



\ + - 1 C + = (2 + 



2 2 
r , m 



/o 2 (l) 



)V'i(l)(l-z m ) + 



2 



3 , /o(l) , ^ 



+ (5 + 



/6'(i) 



2 2 

rim 



)(2 + 



2 2 
rim 



2/6(1) 2/6(l) 2 ' v - ' / 2 (1) 
Solving C + from the above two equations, we have 



/6(i) 2/6(i) 2 
^(i)(i-z m ). 



c 4 



^" A+ (-4/6(1) + m 2 r 2 +2/6(l)z m - 



m 2 r 2 + z m ) 



/6(l)(-2z m - A+ + z m A+) 



(3.21) 



(3.22) 
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Plugging the expression (|3.22| ) into the equation( |3.20|) , we get 



in 



l r 4 (z m -l) fg(l) (m 2 r 



2/^(i) 2 m) 



o ■ -m 2 r 2 + z m \ 



4A + - 2z m A + 

2Zm. ^^4- ~\~ Zm^^^- 



1 1 

m r± 



7 - 6z m + 



[z m - 1)A+ 



+ 



r 2 2 (1 - z m ) 
2/^(l) 2 



- 9 + 7^3.23) 



After plugging = + ^ , = ^ + «Vg and 0' o (l) = - Mo into 



equation Q3.23| ), we obtain 

7-4(162 - 126z m ) m 2 r 4 (36 - 3(k r , 



Z 4 



+m r + (z m - 1) + 



Z 2 

6m 2 r+ A+(l - s m ) 
l 2 ( y —2z m — A+ + z m Aj 



r 4 A + (72 - 36z m ) 
l\-2z m - A + + z m A+) 



+ 



+ re 2 r 2 (54-42z m ) 12/t 2 r 2 A+(z m - 2) { K 2 m 2 r 2 A+(^ m - 1) 



/ 2 



+/t 2 m 2 r 2 (8 - 7z r , 



l 2 ( y —2z m — A + z m A_ 



(7z m -9) | A+(2-z 



2z m — A+ + £ m A+ 

«Vo, (3-24) 



-2z m A-|- ~\~ z m A_j_ 



Neglecting the k 4 terms, we get 



+m 4 / 4 (l - z ; 



(162 - 126z m ) + m 2 / 2 (36 - 30z m ) - 



6m 2 / 2 A+(l - z m ) 
—2z m — A + + z m A_|_ 
12A + (z m -2) 



1 + K 

m 2 l 2 A 



A + (72 - 36^ m ) 
~2z m — A+ + z m A_ 
(54 - 42^ m ) 

I (1 2 m ) 



+ 



l 2 [—2z m — A_|_ + ^ m A_)_)(l — z m ) —2z m — A 4. + z m A_|_ 
m 2 / 2 (8 - 7^ m 



(1 - Zm) 

Further using the relation /i = ■&- , we obtain 



(3.25) 



r+ = v //)/ 2 a/1 - 2n 

+m 4 / 4 (l - z 
4 



(162 - 126z m ) + m 2 / 2 (36 - 30z m ) - 



6m 2 / 2 A+(l -z m ) 
~2z m — A+ + z m A_i_. 
12A + (z m -2) 



1 + k 2 
m 2 Z 2 A 



A + (72 - 36^ m ) 
2z m A_|_ -|- z m A^ 
(54 - 42z m ) 
/ 2 (1 - z m ) 



l 2 {—2z m — A+ + z m A+)(l ~~ ^" 
m 2 Z 2 (8 - 7z m 



-2z m — A+ + z m A_|_ 



(1 ^*m) 

The Hawking temperature in this model is 



T 



3r + 



1 - 



6r 2 



(3.26) 



(3.27) 
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At the critical point, there exists a relationship T = T c . Substituting /xq and r + we 
can obtain 



where 

Air 

i 

T 2 



162 - l2Qz m ) + m 2 Z 2 (36 - 30z„ 
6m 2 / 2 A 11 



(3.28) 



A + (72 - 36z r , 



-2z 



A 



2,21) 



A, 



2/ 2 (l - z m )vl - z r , 



162 - 126^ m ) + m l V (36 - 30z m ) 



+m 4 / 4 (l-0- 



6m 2 / 2 A , (1 - z n 



A+ + 2 m An 



(3.29) 

A + (72 - 36z m ) 
~2z m — A + + £ m A_|_ 

(3.30) 



When setting A + = §, z m = § and m 2 l 2 = -\, we have T c = 0.058^p(l - 5.809^), 
which agrees with the result obtained in ( 2.29Q . Figure § shows that T 2 is positive in 
the range < z m < 1. This means that the effects of the backreaction can make the 
condensation harder to be formed. This result also agrees with the results obtained 

in n m m 0, n . 




Figure 2: T 2 is positive for an arbitrary value of z m . Note that A- 



| and m 2 l 2 



3.2 The upper critical magnetic field with backreaction 

In this part, we study the behavior of a d-wave holographic superconductor which is 
exposed to an external magnetic field and simultaneously consider the backreaction. 



As pointed out in ||T], one may regard the scalar field ip as a perturbation and 
examine its behavior in the neighborhood of the upper critical magnetic field _B c2 . 
In this case, ip is a function of the bulk coordinate z and the boundary coordinates 
(x,y) simultaneously. According to the AdS/CFT correspondence, if the scalar field 
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if) ~ X(x,y)R(z), the vacuum expectation values < O >oc X(x,y)R(z) at the 
asymptotic AdS boundary (i.e. z — » 0) [[33], [43| . One can simply write < O >oc R(z) 
by dropping the overall factor X(x,y). To the leading order, we have 



if)! = ipx(x,y,z), A t = (f) (z), A x = 0, A y = B c2 x. 



(3.31) 



One can check that to the leading order, the ansatz given above is proper for the 
following calculation. Moreover, the black hole carries both electric and magnetic 
charge and the bulk Maxwell field implies that 



A = B c2 xdy + (f) dt. 



(3.32) 



The constant external magnetic field may also backreact on the background black 
hole geometry The tt and rr components of the Einstein's equations in the presence 
of the magnetic field at the leading order become 



w , fo 3r 2 
fo + - ~ -j2 + K r 



-<j)' 2 + m z r + fo[f + 



,2/2 



'.'2 



q 2 <f) 2 ifj 2 e x 



+ 



3t 

4 r 2 



0. 



X' + 2k 2 t if)' 2 + 



q 2 (f) 2 if) 2 e x 
~lo^~ 



Note that near the critical temperature, we have expanded f(z) and x( z ) as 

f = fo + e 2 h + e i U + ... 
X = e 2 X2 + e 4 X4 + ••• 



(3.33) 



(3.34) 
(3.35) 



It is worth noting that the background metric is the dyonic black hole in Ad S4 pi 
because of the magnetic field. To the zeroth order, / is solved as 



fo 



272 



z 2 l 



1 + z + zr 



K 2 l 2 H 2 Z 3 K 2 l 4 B 2 2 Z 3 



2r 2 



2ri 



and the Maxwell fields are given by 

0o = A*o " 



— z. A,, 



B c2 x. 



(3.36) 



(3.37) 



For more than second order 0(e 2 ) , the spacetime metric and the matter fields 
should depend on the spatial coordinates (x, y) and the equations of motion become 
nonlinear partial differential equations. We restrict ourself to the first order and 
indeed this is enough to determine the critical temperature and the upper critical 
magnetic field. Similarly, we have expanded if) and as series in e 



= 0o + e 2 02 + e 4 04 + ••• 
if) = eip! + e 3 ip 3 + e 5 ip 5 + 



(3.38) 
(3.39) 
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At the first order, the equation of motion is 



7o z 



z A f 2 zf 



z 2 fo 



d 2 x + (d y - iB c2 x) 2 Vi- (3.40) 



We assume a separable form for = e lkyV X n ( y x)R n (z), the equation can be 

brought to the form, 



-X'n{x) + [k y - iB c2 x) 2 X n (x) = X n B c2 X n 



x 



(3.41) 



2 2 

r±m 



fo z' 



z A fl z*f ■ zf 



(3.42) 



where A ra stands for the Landau energy level of the harmonic oscillator equation. 
Taking account that the equation ( |3.41| ) is solved by Hermite polynomials and has a 

(k y -iB c2 x) 

series roots like X n (x) = e 2B ^ H n (x), so the solution of ipi can be written as, 



il> 1 = R (z) J %2c j e ik "X (x). 



(3.43) 



We work on the lowest mode n = in what follows, which is the first to condensate 
and is the most stable solution after condensation. Next, we focus on solving the 



equation ( 3.42 ) by the matching method and study the relationship between the 
critical temperature and the magnetic field away from the probe limit. Regularity 
at the horizon requires 



r 2 m 2 _ B c2 l 2 



2+ 



(3.44) 



Near the AdS boundary, _R can be written as Rq(z) = C + z A+ . We expand the 
i?o near the horizon, 



R (z) = R (z)(l) - i2o(*)'(l)(l + -Ro(z)"(l)(l ~ zf + 



(3.45) 



Similarly, we can write the expression of Rq(1) from ( |3.42| ) 

Be/ 2 



K(i) 



2 2 
m r 



2 V 2/5(1) 2/5(1) 2/5(1) 



Rod)' 



2a2 



+ 3 + 



/i(i) m 1 2/5(1)= 



+ 



iUi). 



(3.46) 
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Figure 3: The parameter k 2 shown as a function of T/T c . The plot shows analytic result 
in which Lines from top to bottom are for B C 2 = 100, B C 2 = 1000, B C 2 = 2000. 




Figure 4: The magnetic field B C 2 shown as a function of T/T c . The plot gives analytic 
result for k 2 = (red), k 2 = 0.0001 ( green), k 2 = 0.0004 (blue) from bottom to top. 



Following the above method, we obtain 



= -B 2 c2 + B 2 c2 z m + (-162 - 36m 2 - m 4 + 126z m + 30m 2 z m + m 4 z n 



r 2 (/x 2 - 30B c2 - 2m 2 B c2 - fi 2 z m + 24B c2 z m + 2m 2 B c2 z 
' \ . n .in _ a- i .... i 

?3 

2 1/1-/1 /o. \ d2 



-A_|_ — 2z m + A+z m 



+ k s (54 - 42z m )5 c 2 2 + (8 - 7z m )m 2 B, 



r^.A + i? C 2(6 — 6^ m ) + r 4 _A + (72 + 6m 2 — 36z m — 6m 2 z ( 



1 - z m )A + B 3 c2 



c2 



r 



+ 



-A. 



(24 + m 2 - 12;z m - m 2 z m )A + B 2 c2 2 

h /i 



+ 



— A + — 2z m + A + z m 
I - z m ) A + B c2 (24 + m 2 - 12z m - m 2 z m )r 2 + A + 



(54 + 8m 2 )r 2 + (7 - 6z m )B c2 - (42 + 7m 2 )r 2 z 



-A. 



-A. 



2z m + A + z r , 



(3.47) 
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The difference between (|3.24 ) and ( 3.47 ) comes from the B c2 related terms, from 
which we can derive the relationship between B& and r + as follows 



B 2 c2 



323761r1 16853378855/? 2 r1 35278 



196 



19208 



49 



2 2 2 

k r + fj, Q . 



We find 



T c = 0.058^(1 - 143.463-^) 



(3.48) 



(3.49) 



which means that the critical temperature in external magnetic field is lower than 
the one without magnetic field. We stress that (|3.48|) is not enough to determine 
the relation among the upper critical magnetic field, the system temperature T and 
the critical temperature T c . Of course, there is an ambiguity in the choice of the 
matching radius. But, the result turns out to be fairly insensitive to the choice of 

it. First substitute ( |3.25|) and ( p.48| ) into the Hawking temperature T = ^3 — 



K 2 l 2 fJ? 

2r'i 



2 ,a 



and we see that the Hawking temperature plays the role of the 



critical temperature in the presence of magnetic fields. Taking advantage of the 
relation fM — and equation ( |3.47| ), we finally obtain 



B r 



47r 2 r 6 2 

63 




218617 



105144— 

T-i4 



+ ?CT,T C ), 



(3.50) 



where 



q(T,T c ) = 31231ttV 



31T 2 
441T 2 



105144^ + 6949 



+ 



569T 2 



5292J218617 + 



105144T* 5292T 2 



The result (|3.50| ) implies that it is only applicable near the critical temperature T c 
because the k 2 term will be divergent in the low temperature limit. One may find 
that when k 2 = 0, the result exactly agrees with [fiT] , which is also consistent with 

the Ginzburg-Landau theory where B c2 oc ^1 — T/T^j . Figure § shows that for fixed 

B c2 , the critical temperature drops as k 2 increases. As shown in Figure (|, the upper 
critical magnetic field decreases as T/T c goes up and vanishes at T = T c . We also 
find that the coefficient of the k 2 term is positive for the system temperature T (see 
Fig. H). This result also agrees with the result obtained in 



4. Numerical results 



In this section, we will visually show the relationship among the critical temperature, 
the upper critical magnetic field and the backreaction parameter k 2 by using the 
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0.5 0.6 0.7 0.8 0.9 1.0 

T 



Figure 5: The coefficient of the k 2 term of the upper critical magnetic field as a function 
of the temperature T/T c . 

numerical shooting method. We set r + = 1 and I = 1 in the numerical computation. 
Firstly, we want to show how the critical temperature varies as k 2 changes. From 
Figure [], one can see that for fixed the critical temperature decreases as k 2 
rises up. That is to say, when applying the external magnetic field, the backreaction 
of the spacetime makes condensation harder to happen. This is consistent with the 
analytic result showed in Figure |3|. 




0.020 0.025 0.030 0.035 

T 

T c 



Figure 6: The numerical result in which lines from top to bottom are for B c2 = 100 (red), 
B c2 = 200 (green), B c2 = 300 (blue). 

Secondly, we plot the curves between the critical temperature and the upper 
critical magnetic field B c2 for fixed parameter k 2 . When setting different values 
of backreaction n 2 = 0, n 2 = 0.01, k 2 = 0.04 in Figure 0. We find the magnetic 
field B c2 goes to the opposite direction with increasing T/T c . The numerical result 
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also demonstrates that when k 2 = 0, the upper critical magnetic field B c2 is finite 
and smaller than k 2 ^ cases. That is to say, although the critical temperature is 
significantly suppressed by a non-zero k 2 , the upper bound of B c2 become larger. 

The numerical computation is done by using the shooting method. The plot of 
the analytic result (Figures |3] and |j) depends on the choice of the matching point. 
But the analytic results given in Figure |3] and f| are qualitatively in good agreement 
with the numerical results shown in Figures ^| and [7]. 

200 
150 
B c2 100 
50 


0.90 0.92 0.94 0.96 0.98 1.00 

T 

T c 

Figure 7: The numerical result for Lines from top to bottom are for k 2 = (red), k 2 = 0.01 
(green), k 2 = 0.04 (blue) respectively. 

5. Conclusions 

In this work we studied the d-wave holographic superconductors in the probe limit 
and away from the probe limit by using the analytic and numerical computation. 
In the probe limit, we obtained analytic expressions for the order parameter, the 
critical temperature and the upper critical magnetic field. The analytic calculation 
is useful for gaining insight into the strong interacting system. The result shows that 
the upper critical field gradually increases to its maximum value B c2 at absolute 
zero temperature T = 0, while vanishing at the critical temperature T = T c . This 
behavior satisfies the relation given in the Ginzburg-Landau theory. Away from 
the probe limit, we found that the backreaction of the electric field results in the 
AdS-Reissner-Nordstrom black hole solution, which in turn hinder the formation of 
the condensation. We obtained the analytic formula for the critical temperature 
in this case and found T 2 is always positive. Condensation becomes even harder 
to happen in the presence of a constant external magnetic field. The background 
spacetime geometry to the leading order becomes a dyonic black hole solution in AdS 
space. The coefficient of spacetime backreaction on the upper critical magnetic field 
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is positive, which indicates that the magnetic field becomes strong with respect to 
the backreaction and is in consistent with the result observed in |p4| |. We have also 
shown the corresponding numerical results for each case. Qualitatively our analytic 
and numerical computational results matched very well. 

The method employed in this work is very similar to that used in jTI]. But 
surprisingly, we found that the spacetime backreaction presents us an interesting 
property of holographic superconductors: While the backreaction causes the sup- 
pression of the critical temperature, it can enhance the upper critical magnetic field 
because the 0(k 2 ) term is positive. Moreover, we have checked the analytic result by 
comparing with numerical computation. In fact, we found that the holographic su- 
perconductors with backreaction is very similar to carbon doping in MgB2 reported 
in recent experiments ||45|| : where with the results of the suppression in T c , while the 
B C 2 performance is improved. The consistency of the phenomenon in the external 
magnetic field in holographic superconductor and MgB 2 superconductor is interesting 
and further understanding on their relation is called for. 
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